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A  theoretical  model  of  solid  oxide  electrolysis  cells  considering  the  heterogeneous  elementary  reactions, 
electrochemical  reactions  and  the  transport  process  of  mass  and  charge  is  applied  to  study  the  relative 
performance  of  H20  electrolysis,  C02  electrolysis  and  C02/H20  co-electrolysis  and  the  competitive 
behavior  of  heterogeneous  chemical  and  electrochemical  reactions.  In  cathode,  heterogeneous  chemical 
reactions  exist  near  the  outside  surface  and  the  electrochemical  reactions  occur  near  the  electrolyte. 
According  to  the  mathematical  analysis,  the  mass  transfer  flux  DVc  determines  the  main  zone  size  of 
heterogeneous  chemical  reactions,  while  the  charge  transfer  flux  aWV  determines  the  other  one.  When 
the  zone  size  of  heterogeneous  chemistry  is  enlarged,  more  C02  could  react  through  heterogeneous 
chemical  pathway,  and  polarization  curves  of  C02/H20  co-electrolysis  could  be  prone  to  H20  electrolysis. 
Meanwhile,  when  the  zone  size  of  electrochemistry  is  enlarged,  more  C02  could  react  through  elec¬ 
trochemical  pathway,  and  polarization  curves  of  C02/H20  co-electrolysis  could  be  prone  to  C02  elec¬ 
trolysis.  The  relative  polarization  curves,  the  ratio  of  C02  participating  in  electrolysis  and  heterogeneous 
chemical  reactions,  the  mass  and  charge  transfer  flux  and  heterogeneous  chemical/electrochemical  re¬ 
action  main  zones  are  simulated  to  study  the  effects  of  cathode  material  characteristics  (porosity,  particle 
diameter  and  ionic  conductivity)  and  operating  conditions  (gas  composition  and  temperature). 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

CO2/H2O  co-electrolysis  in  solid  oxide  electrolysis  cells 
(SOECs)  is  one  of  the  efficient  routes  to  reduce  CO2  emission 
and  to  convert  and  store  the  renewable  power  into  the  form  of 
the  chemical  energy  of  hydrocarbon  fuel  1-3  .  In  the  cathode, 
H2O  and  CO2  electrochemical  reductions  occur  at  the  triple 
phase  boundary  (TPB).  The  reversible  water  gas  shift  (WGS  and 
RWGS)  reactions  also  react  on  the  electrode  catalyst  (such  as 
nickel): 
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H20  +  2e~  ->  H2+02_ 

(i) 

C02+2e~  -»CO  +  02~ 

(2) 

RWGS 

h2+co2^^h2o  +  co 

WGS 

(3) 

It  is  generally  believed  that  CO  can  be  reduced  from  CO2  elec¬ 
trolysis  and  from  RWGS  reaction,  but  it  is  not  clear  that  how  much 
CO  is  produced  from  RWGS  reaction  [3-5  .  Two  different  conclu¬ 
sions  can  be  found  in  the  published  experimental  studies  of  CO2 / 
H2O  co-electrolysis:  (1)  CO2/H2O  co-electrolysis  performs  almost 
the  same  as  H2O  electrolysis,  CO  is  nearly  all  generated  from  RWGS 
reaction  [6,7];  (2)  C02/H20  co-electrolysis  performs  between  H20 
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electrolysis  and  CO2  electrolysis,  CO  is  partial  generated  from  RWGS 
reaction  [3,8  .  Ni  [2,9  has  developed  both  one-dimensional  and 
two-dimensional  thermal-fluid  model  and  studied  that  the  per¬ 
centage  of  CO  generated  from  RWGS  reaction  highly  depends  on 
the  operating  temperature  and  inlet  gas  composition  of  cathode. 

Except  for  temperature  and  gas  composition,  we  believe  the 
mass  and  charge  transport  process  in  cathode  also  play  an  impor¬ 
tant  role  on  the  CO2  reaction  pathway  and  the  relative  performance 
of  H2O  electrolysis,  CO2  electrolysis  and  CO2/H2O  co-electrolysis. 
The  cathode  thickness  which  influences  the  transport  length  has 
been  proved  to  be  one  of  the  factors,  through  establishing  and 
validating  an  elementary  reaction  model  in  our  previous  work  10]. 
It  is  found  that  the  rate  of  heterogeneous  chemical  reactions  is 
about  20-100  times  faster  than  that  of  electrochemical  reactions, 
and  the  two  kinds  of  reactions  occur  in  the  different  zones  in 
cathode  as  shown  in  Fig.  1 :  the  heterogeneous  chemical  reactions 
mainly  react  near  the  outside  surface  of  cathode,  while  the  elec¬ 
trochemical  reactions  mainly  react  near  the  electrolyte.  It  could  be 
speculated  that  at  the  same  local  reaction  site,  if  the  reactant  gas 
doesn't  reach  the  equilibrium,  the  gas  prefers  to  heterogeneous 
chemical  reactions  rather  than  electrochemical  reactions.  Thus,  for 
a  thick  cathode,  the  electrochemical  reaction  main  zone  is  far  from 
that  of  heterogeneous  chemical  reaction,  so  that  there  are  sufficient 
reaction  sites  for  gas  reacting  first  through  chemical  pathway  to 
reach  the  equilibrium  and  then  through  electrochemical  pathway. 
For  CO2/H2O  co-electrolysis,  CO2  can  both  react  through  chemical 
and  electrochemical  reaction  pathways,  and  the  electrochemical 
performance  lies  between  H2O  and  CO2  electrolysis.  If  the  cathode 
is  relatively  thin,  the  two  zones  will  overlap  and  compete  with  each 
other.  Finally,  the  gas  prefers  to  participate  in  chemical  reactions 
first  and  most  of  CO2  reacts  through  chemical  reaction  pathway, 
which  causes  the  performance  approaching  to  H2O  electrolysis. 

In  this  paper,  the  objective  is  to  further  study  the  impacts  of  mass 
and  charge  transport  on  the  relative  performance  of  H2O  electrolysis, 
C02  electrolysis  and  C02/H20  co-electrolysis  and  the  competitive 
behavior  of  heterogeneous  chemical  and  electrochemical  reactions, 
which  will  be  beneficial  for  deeper  understanding  of  the  cathode 
mechanism  of  CO2/H2O  co-electrolysis,  cell  design  and  optimizations. 


Table  1 

Heterogeneous  reaction  mechanism  on  Ni-based  catalysts. 


Reaction 


A  (cm,  mol,  s)a  na  E 

(kj  mol_1)a 


Adsorption  &  desorption 

1 f  H2  (g)  +  (Ni)  +  (Ni)  -  H(Ni)  +  H(Ni) 
lr  H(Ni)  +  H(Ni)-(Ni)  +  (Ni)+H2(g) 
2f  02(g)  +  (Ni)  +  (Ni)-O(Ni)  +  O(Ni) 
2r  O(Ni)  +  O(Ni)  -  (Ni)  +  (Ni)+02  (g) 
3f  H20(g)  +  (Ni)  — ►H20(Ni) 

3r  H20(Ni)  — >  (Ni)+H20(g) 

4f  C02(g)  +  (Ni)  — >C02(Ni) 

41  C02  (Ni)  — ►  (Ni)+C02  (g) 

5f  CO(g)  +  (Ni)->CO(Ni) 

51  CO(Ni)  —>  (Ni)  +  CO(g) 

Surface  reactions 

6f  H(Ni)  +  O(Ni)  -+OH(Ni)  +  (Ni) 

61  OH(Ni)  +  (Ni)  — >H(Ni)  +  O(Ni) 

7f  H(Ni)  +  0H(Ni)-»H20(Ni)  +  (Ni) 

T  H20(Ni)  +  (Ni)^H(Ni)  +  OH(Ni) 

8f  OH(Ni) +  0H(Ni)^H20(Ni) +  0(Ni) 
8r  H20(Ni)  +  O(Ni)  — >OH(Ni)  +  OH(Ni) 
9f  CO(Ni)  +  0(Ni)^C02(Ni)  +  (Ni) 

91  C02  (Ni)  +  (Ni)  CO(Ni)  +  O(Ni) 


1.000 

X 

10  02b 

0.0 

0.00 

2.545 

X 

10+19 

0.0 

81.21 

1.000 

X 

10~°2b 

0.0 

0.00 

4.283 

X 

10+23 

0.0 

474.95 

0.100 

X 

10  0°b 

0.0 

0.00 

3.732 

X 

10+12 

0.0 

60.79 

1.000 

X 

10  05b 

0.0 

0.00 

6.447 

X 

10+07 

0.0 

25.98 

5.000 

X 

10~°2b 

0.0 

0.00 

3.563 

X 

10+11 

0.0 

111.27 

^CO(s) 

-50.00 

5.000 

X 

10+22 

0.0 

97.90 

1.781 

X 

10+21 

0.0 

36.09 

3.000 

X 

10+2° 

0.0 

42.70 

2.271 

X 

10+21 

0.0 

91.76 

3.000 

X 

10+21 

0.0 

100.00 

6.373 

X 

10+23 

0.0 

210.86 

2.000 

X 

10+19 

0.0 

123.60 

^CO(s) 

10+23 

-50.00 

4.653 

X 

-1.0 

89.32 

a  Arrhenius  parameters  for  the  rate  constant  written  in  the  form:  k  =  ATnexp(-£/ 
RT). 

b  Sticking  coefficient. 
c  Coverage  dependent  activation  energy. 


The  elementary  reaction  model  of  SOECs  developed  and  validated  in 
our  previous  work  [10]  continues  to  be  used  in  this  simulation.  The 
sizes  of  the  heterogeneous  chemical/electrochemical  reaction  main 
zones  are  quantified  according  to  the  mathematical  analysis.  The 
effects  of  material  characteristics  of  cathode  and  operating  conditions 
of  SOECs  are  studied  in  details.  Here,  the  microstructure  including 
porosity  cCa  and  particle  diameter  dca,  and  the  ionic  conductivity 
^ion.ca  are  chosen  as  the  typical  parameters  of  material  characteristics. 
In  addition,  as  the  important  operating  parameters,  the  effects  of 
temperature  T  and  gas  composition  are  also  discussed. 
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Fig.  1.  Heterogeneous  chemical/electrochemical  reaction  main  zones,  C02  reaction  pathway  and  mathematical  analysis  of  the  two  zones. 
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Table  2 

Summary  of  the  governing  equations  of  the  model.  Details  can  be  found  in  Refs.  [10—12]. 


Cathode  heterogeneous  chemistry 

The  reaction  rate  constant  of  surface  reactions  and 
desorption  reactions  the  Arrhenius  form 

The  reaction  rate  constant  of  adsorption  reactions 


Cathode  electrochemistry 

The  current  source 

The  forward  electrochemical  reaction  rates 


The  reverse  electrochemical  reaction  rates 


The  cathode  overpotential 

Reaction  area 

The  TPB  active  area  per  unit  volume 


The  effective  Ni  surface  area  per  unit  volume 


The  total  number  of  particles  per  unit  volume 

Charge  balance 

The  ionic  charge  equation  in  cathode 
The  electronic  charge  equation  in  cathode 
The  ionic  charge  equation  in  anode 


The  electronic  charge  equation  in  anode 
The  ionic  charge  equation  in  electrolyte 

The  anode  overpotential 

Mass  balance 

The  mass  balance  equation  in  porous  electrode 
The  effective  diffusivity  of  gas  species  k 


The  effective  molecular  diffusion  coefficient 


The  effective  binary  molecular  diffusion  coefficient 


The  effective  Knudsen  diffusion  coefficient 


The  source  terms  of  both  gas  and  surface  species  in  fuel  electrode 


ki  =A  Jn'  exp(  -JM  Yl  eT  exP 


I<g+Ks 


RT, 


k=  1 


£ki  f_k 

RT 


ki  =  aiTb,ex p(-A)  T 


/  RT 
2tzW 


(4) 


(5) 


Qca  =  2F^/<ecC0(Ni)c(YSZ)  “  ^~ecC02- (YSZ)C(^)  J^IPB 
*0 


/<ec  — 


^TPBco(Ni)C(YSZ) 


exp 


-2(1 -a 


FVca 


RT 


lo 


f5TpBC02-(YSZ)l'(Ni) 

Vca  =  ^elec,ca  _  ^ion,ca  _  ^ref. 


ca 


(6) 

(7) 

(8) 

(9) 


Trsin2  dr^pntnepnipZepZipPepPip 


^TPB  = 


5Ni  =  TCrepntnep  (  4  - 


z 

sin  0nipZepZip  sin  ^rienZenZe 


lep^ep^ep 


(10) 


(ID 


nt  = 


—  £ 


3TCrep(nep  +n[pa^) 


(12) 


(  °fon,caV^ion,ca  )  —  Qion,ca 


V-  - 


^elec.ca^^elec^ca  )  —  Qelec,ca  —  Qion.ca 


(13) 


(14) 


V-  - 


6ff  T7\/ 
°ion,an  v  v ion, an 


=0, 


(3RT 


on.an 


exp  - 


exp 


4  F 

2(1  -  a)F7]an 


‘-air 

RT 


air 


0.25 


5tpb 


an 


RT 


^  ’  (  °elec:an^^elec,an  )  —  Qelec,an  —  Qion,an 


^an  —  ^elec,an  ^ion,an  ^ref,an 


Cof  c_p^an' 
-bulkeXP(  RT 

'O,  v 


(15) 

(16) 

(17) 

(18) 


V(  —  Df/fVC/,  rj  )  =R 


Df 


k  VL/cg)=^/cg 


+ 


-1 


Deff 

uk.  mole 


Deff 
uk.  Kn 


(19) 


(20) 


D 


,eff 

/c.mole 


1  -xk 


I) 


(21) 


Deff  0.0010UT1  75(1/Mt  +  1/M,)1/2 
,<J  rp(v^  +  v^)2 

neff  [8RT_ 

k’Kn  3  T  V  ^Mk 

N+2  /  I<g+I<s+2 

i=l  k= 1 

D  Qelec,an 


(22) 


(23) 

(24) 

(25) 


The  source  terms  in  anode 
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Fig.  2.  Simplified  modeling  geometry  and  materials. 


2.  Model  development,  experiment  and  model  validation 

A  one-dimensional  elementary  reaction  model  of  SOECs  using 
the  finite  element  commercial  software  COMSOL  Multiphysics  is 
developed  and  described  in  details  in  our  previous  work  10-13]. 
This  model  couples  with  elementary  heterogeneous  reactions 
[14,15  ,  electrochemical  reactions  16,17  ,  charge  balance  [18-20] 
and  mass  balance  considering  molecular  diffusion  and  Knudsen 
diffusion  14,21-23].  A  heterogeneous  reaction  mechanism  on  Ni- 
based  catalysts  simplified  from  the  works  of  Deutschmann’s 
group  [15,24 1  is  used  and  listed  in  Table  1.  The  effective  TPB  area 
and  Ni  surface  area  are  calculated  by  the  particle  coordination 
number  theory  in  binary  random  packing  of  spheres  and  the 
percolation  theory  20,21].  The  governing  equations  of  the  model 
are  summarized  in  Table  2. 

Cathode  supported  button  cells  made  by  Shanghai  Institute  of 
Ceramics  Chinese  Academy  of  Sciences  were  tested  by  a  button  cell 
reactor  and  measurement  system  [11,12].  The  simplified  modeling 
geometry  and  materials  are  shown  in  Fig.  2.  This  SOEC  model  was 
validated  with  the  polarization  curves  of  H2O  electrolysis,  CO2 
electrolysis  and  CO2/H2O  co-electrolysis  with  different  composi¬ 
tions  at  700  °C  10].  In  this  paper,  all  the  model  parameters  are  the 
same  as  the  Ref.  [10]  and  would  not  be  changed. 

3.  Simulation  conditions  and  parameters  definition 

Table  3  summaries  the  main  simulation  conditions  used  in  the 
model.  The  cathode  porosity  eca  and  cathode  particle  diameter  dca 


Table  3 

Simulation  conditions  of  the  model. 


Parameter 

Value,  Unit 

Microstructure 

Porosity  of  cathode, 

eca 

Experiment  value 

0.335  [10] 

Simulation  value 

0.3— 0.6 

Particle  diameter  of  cathode,  dca  (assumption:  particle  diameter  = 

pore  diameter) 
Experiment  value 

0.60  [10]  pm 

Simulation  value 

0.6— 1.2  pm 

Ionic  conductivity  of  cathode,  <Ti0n,ca 

Validation  value 

YSZ:  3.34  x  104- 

exp  (— 10300/r1)  [25]  S  itT1 

ScSZ:  6.92  x  104 

•exp  (-9681/7  [11]  S  m1 

Simulation  value 

Gas  composition 

Validation  value 

x  1. 0-4.0 

Cathode 

Oxidant:  reductant  H20/H2 

h2o/h2/co2/co 

co2/co 

4:1 

80/20 

40/10/40/10 

80/20 

3:2 

60/40 

30/20/30/20 

60/40 

1:1 

50/50 

25/25/25/25 

50/50 

2:3 

40/60 

20/30/20/30 

40/60 

1:4 

20/80 

10/40/40/10 

20/80 

Anode 

02/N2  =  21/79 

Temperature 

700-900  °C 

Pressure 

1.0  bar 

are  changed  in  the  ranges  of  0.3-0.6  and  0.6-1.2  pm,  respectively. 
Here,  the  particle  diameters  of  ionic  and  electronic  conductors  in 
electrodes  are  assumed  to  be  the  same  and  equal  to  the  pore 
diameter  23,26  .  The  cathode  ionic  conductivity  <7ion,ca  is  artificially 
multiplied  by  a  factor  ranged  from  1.0  to  4.0.  For  comparison,  the 
cathode  gas  compositions  of  H2O  electrolysis,  CO2  electrolysis  and 
CO2/H2O  co-electrolysis  are  given  as  the  same  ratio  of  oxidant  and 
reductant  which  is  varied  from  4:1  to  1 :4.  The  ratio  of  H2O  and  CO2 
and  ratio  of  H2  and  CO  for  co-electrolysis  are  all  fixed  as  1.  Air  is 
supplied  in  the  anode.  The  temperature  is  changed  from  700  °C  to 
900  °C,  while  pressure  is  kept  at  1.0  bar. 

Quantifying  the  size  of  the  heterogeneous  chemical/electro¬ 
chemical  reaction  main  zones  could  supply  a  more  legible  and  vi¬ 
sual  way  to  analyze.  As  shown  in  Fig.  1  and  Table  2,  the 
heterogeneous  chemistry  is  mainly  calculated  by  the  mass  balance 
equation  (Eq.  (19)),  and  the  information  of  reaction  rates  are  all 
inside  the  sources  R  (Eq.  (24)).  Thus,  the  variation  speed  of  c,  the 
mass  transfer  flux  DVc,  could  signify  the  overall  reaction  rate  of 
heterogeneous  chemistry.  Similarity,  the  electrochemistry  is  mainly 
calculated  by  the  charge  balance  equation  (Eq.  (13)),  and  its  source 
Q  involves  in  the  electrochemical  reaction  rates  (Eqs.  (6)-(8)).  So 
the  charge  transfer  flux  a  W,  the  term  in  the  parentheses  of  Eq.  (13), 
denotes  the  overall  reaction  rate  of  electrochemistry.  We  artificially 
define  that  the  boundary  of  main  zone  of  chemical  reactions  is  the 
location  where  c  has  been  90%  changed  without  current  (at  OCV), 
meanwhile  the  boundary  of  electrochemical  reactions  is  the  loca¬ 
tion  where  V  has  been  90%  changed  with  current  (for  example  at 
1.4  V).  This  specific  value  90%  doesn't  have  any  physical  meaning 
here,  and  can  be  any  other  values  like  95%  or  85%.  This  value  is  only 
used  to  quantify  the  zone  size  and  further  to  compare  and  analysis 
its  variation. 

The  mass  transfer  flux  for  all  the  gas  species  H2O,  H2,  CO2,  CO  is 
all  the  same  on  the  basis  of  mathematics.  Besides,  according  to  the 
gas  distribution  in  cathode  calculated  in  Fig.  6  of  Ref.  10],  it  was 
shown  that  the  heterogeneous  chemical  reactions  on  Ni  product 
CO2  and  H2  but  consume  CO  and  H2O  at  700  °C;  the  electrochemical 
reactions  product  CO  and  H2  but  consume  CO2  and  H2O.  Thus,  CO2 
and  CO  show  the  opposite  behaviors  on  these  two  kinds  of  re¬ 
actions,  and  the  concentration  of  CO2  is  chosen  to  quantify  the  zone 
here.  As  shown  in  Fig.  3(a),  the  heterogeneous  chemical  reaction 
main  zone  is  quantified  as  the  zone  from  the  outside  surface  of 
cathode  to  the  location  with  90%  AcC02  at  OCV.  On  the  other  hand, 
the  ionic  transport  is  much  slower  than  electronic  transport  and 
limits  the  charge  transport,  so  that  the  ionic  potential  is  used  to 
quantify  the  zone.  The  electrochemical  reaction  main  zone  includes 
the  zone  from  the  electrolyte  to  the  location  with  90%  AVion  at  1.4  V, 
as  shown  in  Fig.  3(b). 

In  order  to  clear  how  much  CO2  reacts  through  chemical  or 
electrochemical  reaction  pathway,  the  ratio  of  the  variation  of  CO2 
concentration  reduced  by  electrolysis  AcC02i  electrochemical  and  the 
variation  caused  by  chemical  reactions  AcC02i  chemical  is  quantified 
as  follows  and  expressed  in  Fig.  3(a): 
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Fig.  3.  Quantification  of  the  sizes  of  heterogeneous  chemical/electrochemical  reaction  main  zones. 


_  ^C02,  electrochemical  _  ^C02,  el  (OCV)  Cc02,  elO-^V) 

^  ACco2,  chemical  ^C02,  el  (OCV)  —  Cqq2  cc 

(26) 

Acco2,  electrochemical  is  approximately  considered  as  the  difference  of 
concentration  at  the  boundary  of  electrolyte  cC02  ei  at  OCV  and 
1.4  V.  Besides,  AcC02i  chemical  could  be  equal  to  the  difference  of 
concentration  at  the  boundary  of  electrolyte  cC02iei  and  concen¬ 
tration  at  the  boundary  of  cathode  surface  cC02  c c  at  OCV. 


4.  Results  and  discussion 

4.2.  Effect  of  microstructure:  porosity  and  particle  diameter 

4.1 .2.  Polarization  curves  and  rer/cr 

Fig.  4  represents  and  compares  the  polarization  curves  of  H2O 
electrolysis,  CO2  electrolysis  and  CO2/H2O  co-electrolysis  with 
different  cathode  porosities  eca  and  particle  diameters  dca.  The 
temperature  is  fixed  as  700  °C,  and  the  ratio  of  oxidant  and 
reductant  is  1:1.  The  electrochemical  performances  of  the  three 


Fig.  4.  Polarization  curves  comparison  of  H20  electrolysis,  C02  electrolysis  and  C02/H20  co-electrolysis  with  different  cathode  porosities  (a-c)  and  different  cathode  particle  di¬ 
ameters  (d-f)  at  700  °C. 
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Table  4 

Ratio  (%)  of  C02  participating  in  electrolysis  and  heterogeneous  reactions  at  1.4  V. 


Porosity, 

eca 

Particle  diameter,  dca 

Ionic  conductivity,  (7ion  cac 

0.30 

5.24% 

0.6  pm 

4.43% 

xl.O 

2.29% 

0.35 

4.14% 

0.7  pm 

3.73% 

xl.5 

2.40% 

0.40 

3.35% 

0.8  pm 

3.24% 

x2.0 

2.47% 

0.45 

2.76% 

0.9  pm 

2.88% 

x2.5 

2.52% 

0.50 

2.29% 

1.0  pm 

2.60% 

x3.0 

2.55% 

0.55 

1.92% 

1.1  pm 

2.39% 

x3.5 

2.58% 

0.60 

1.61% 

1.2  pm 

2.22% 

x4.0 

2.60% 

Oxidant:  reductant 

Temperature,  T6 

4:1 

0.35% 

700  °C 

2.29% 

3:2 

0.74% 

750  °C 

3.60% 

1:1 

2.29% 

800  °C 

5.06% 

2:3 

6.68% 

850  °C 

6.60% 

1:4 

12.03% 

900  °C 

8.23% 

a  dca  =  0.6  pm,  oxidant:  reductant  =  1:1  ,T=  700  °C. 
b  £ca  =  0.335,  oxidant: reductant  =  1:1,  T=  700  °C. 
c  £ca  =  0.5,  dca  =  0.6  pm,  oxidant:reductant  =  1:1,7=  700  °C. 
d  £ca  =  0.5,  dca  =  0.6  pm,  7  =  700  °C. 
e  £ca  =  0.5,  dca  =  0.6  pm,  oxidant:reductant  =  1:1. 

electrolysis  modes  all  decrease  with  the  enhanced  eca  in  the 
range  of  0.30-0.60,  when  dca  is  fixed  as  0.6  pm,  and  also  decrease 
with  the  increasing  dca  from  0.6  to  1.2  pm  at  the  same  eca  0.335. 
According  to  Eqs.  (10)-(12)  in  Table  2,  the  TPB  and  Ni  surface 


areas  are  both  proportional  to  l/dca  and  (1  -  £Ca),so  that  the 
decreasing  variation  is  mainly  attributed  to  the  decreasing  reac¬ 
tion  active  sites  per  unit  volume  when  eca  or  dca  increases. 
However,  to  reduce  porosity  and  particle  diameter  blindly  cannot 
always  improve  performance  due  to  the  gas  transport  limitation 
[27]. 

It  is  more  interesting  to  observe  that  the  microstructure  can 
not  only  affect  the  absolute  cell  performance,  but  also  can  signif¬ 
icantly  affect  the  relative  cell  performance  of  the  three  electrolysis 
modes.  The  results  in  Fig.  4(a-c)  show  that  as  eca  grows,  the  po¬ 
larization  curve  of  CO2/H2O  co-electrolysis  is  gradually  close  to 
that  of  H20  electrolysis.  As  shown  in  Fig.  4(d-f),  increasing  dca  also 
represents  a  similar  tendency  of  the  relative  performance,  though 
the  absolute  performances  all  get  worse.  For  CO2/H2O  co¬ 
electrolysis,  if  only  H2O  would  participate  in  electrolysis  and  CO2 
would  be  inert  in  electrolysis,  the  polarization  curves  would  be 
similar  to  that  for  H2O  electrolysis,  like  the  experimental  results  in 
Refs.  [6,7  .  If  co-electrolysis  performs  between  that  of  H20  and  C02 
electrolysis  at  the  same  ratio  of  oxidant  and  reductant,  it  suggests 
that  H2O  and  CO2  electrolysis  both  occur  in  the  electrode  [3,4]. 
Thus,  we  speculate  that  at  the  same  oxidant  and  reductant  ratio, 
the  tendency  of  co-electrolysis  closing  to  H2O  electrolysis  denotes 
more  H2O  and  less  CO2  could  be  reduced  by  electrolysis  and  vice 
versa. 


0  100  200  300  400  500  600  695 


0  100  200  300  400  500  600  695 

Distance  from  cathode  outside  surface(pm) 


Fig.  5.  Surface  species  concentrations  of  O(Ni)  and  (Ni)  of  H20  electrolysis,  C02  electrolysis  and  C02/H20  co-electrolysis  with  different  cathode  porosities  (a,  b)  and  different  cathode 
particle  diameters  (c,  d)  (1.4  V,  700  °C). 
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Fig.  6.  Effects  of  porosity  (a-c)  and  particle  diameter  (d-f)  on  heterogeneous  chemical/electrochemical  reaction  main  zones  in  cathode  for  C02/H20  co-electrolysis. 


Fig.  7.  Polarization  curves  comparison  of  H20  electrolysis,  C02  electrolysis  and  C02/H20  co-electrolysis  with  cathode  ionic  conductivity  multiplied  by  a  coefficient  ranged  from  1.0  to 
4.0  at  700  °C. 
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The  speculation  could  be  directly  proved  by  the  ratio  of  CO2 
participating  in  electrolysis  and  heterogeneous  reactions  rer/cn  as 
listed  in  Table  4.  As  mentioned  before,  the  reaction  rate  of  hetero¬ 
geneous  chemistry  is  highly  faster  than  electrochemistry,  so  the 
ratio  here  is  a  little  small.  It  can  be  observed  that  rer/c r  decreases 
indeed  with  the  enhanced  porosity  and  particle  diameter,  which 
accords  with  the  tendency  of  relative  polarization  curves.  rer/c r  at 
eca  =  0.60  is  only  about  30%  of  that  at  eca  =  0.30,  while  rer/c r  at 
dca  =  1.2  pm  is  about  50%  of  that  at  dca  =  0.6  pm.  Therefore,  these 
results  indicate  that  when  increasing  porosity  or  particle  diameter, 
the  electrochemical  reaction  of  CO2/H2O  co-electrolysis  is  prone  to 
H20  electrolysis,  and  more  C02  reacts  through  the  heterogeneous 
chemical  reaction  pathway. 

4.1.2.  O(Ni)  and  (Ni)  concentrations 

According  to  our  previous  work  TO],  as  the  main  surface  species 
of  the  elementary  electrochemical  reaction  in  this  model  shown  in 
Reaction  (27),  the  surface  species  concentrations  of  O(Ni)  and  (Ni) 
are  proven  to  be  one  of  the  influencing  factors  which  control  the 
relative  performance  of  the  three  electrolysis  modes,  after  studying 
all  the  rate  of  electrochemical  reaction  and  heterogeneous 
elementary  reactions  in  details.  O(Ni)  denotes  the  O  element 
covered  on  the  Ni  surface,  and  (Ni)  denotes  the  free  surface  active 
sites  on  the  Ni  surface. 


zone  is  enlarged  from  110  pm  to  178  pm,  while  the  size  of  elec¬ 
trochemical  zone  slightly  increases  only  3  pm. 

From  Fig.  6(a)  and  (c),  the  mass  transfer  flux  D^2VcCq2  falls 
down  quickly  with  the  distance  from  outside  surface  of  cathode  and 


O(Ni)  +  (YSZ)+2e“  ;  =  02“  (YSZ)  +  (Ni)  (27) 

k~ec 

Fig.  5  shows  the  surface  concentrations  of  O(Ni)  and  (Ni)  at  a 
certain  cell  voltage  of  1.4  V.  In  the  zone  near  the  outside  surface  of 
cathode,  the  variations  of  both  O(Ni)  and  (Ni)  for  CO2/H2O  co¬ 
electrolysis  are  quite  bigger  than  the  other  two  modes  due  to  the 
heterogeneous  chemical  reactions.  Moreover,  the  porosity  and 
particle  diameter  also  both  have  bigger  influence  on  the  O(Ni)  and 
(Ni)  concentrations  in  co-electrolysis  mode  than  that  in  the  other 
two  modes.  When  eca  or  dca  increases,  the  average  O(Ni)  concen¬ 
tration  of  CO2/H2O  co-electrolysis  increases  and  the  average  (Ni) 
concentration  decreases,  and  both  of  their  curves  become  close  to 
the  curves  of  H2O  electrolysis.  Thus,  although  the  absolute  forward 
and  reverse  electrochemical  reaction  rates  (ke c,  k-ec)  of  the  three 
modes  all  decrease,  the  relative  electrochemical  performance  of 
CO2/FI2O  co-electrolysis  is  prone  to  that  of  H2O  electrolysis,  after 
multiplying  the  rate  by  the  surface  concentrations  of  O(Ni)  and  (Ni) 
(Eq.  (6)).  ' 


4.1.3.  The  sizes  of  heterogeneous  chemical/electrochemical  reaction 
main  zones 

O(Ni)  and  (Ni)  concentrations  offer  an  explanation  at  a  micro 
and  elementary  level,  moreover,  the  analysis  of  heterogeneous 
chemical/electrochemical  reaction  main  zones  can  interpret  this 
phenomena  at  a  macro  and  more  legible  way. 

Fig.  6  represents  the  effects  of  porosity  and  particle  diameter  on 
mass  and  charge  transfer  flux  which  are  the  determinant  param¬ 
eters  of  the  size  of  two  zones  as  analyzed  in  Section  3,  and  sum¬ 
marizes  the  sizes  of  heterogeneous  chemical/electrochemical 
reaction  main  zones  in  a  bar  chart.  The  blue  (in  the  web  version) 
column  denotes  the  heterogeneous  chemical  reaction  main  zone 
while  the  green  (in  the  web  version)  denotes  the  electrochemical 
one.  The  results  indicate  that  the  microstructure  strongly  affects 
the  main  zone  size  of  heterogeneous  chemical  reaction,  but  slightly 
affects  that  of  electrochemical  reaction.  The  size  of  chemical  zone 
increases  from  101  pm  to  188  pm  when  eca  increases  from  0.30  to 
0.60,  while  the  size  of  electrochemical  zone  is  only  reduced  20  pm. 
When  dca  is  varied  in  the  range  of  0.6-1.2  pm,  the  size  of  chemical 


■  Heterogenous  chemical  ■  Electrohemical 

reaction  main  zone  reaction  main  zone 
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Fig.  8.  Effects  of  cathode  ionic  conductivity  on  heterogeneous  chemical/electro- 
chemical  reaction  main  zones  in  cathode  for  C02/H20  co-electrolysis. 
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then  gradually  approaches  to  zero,  without  the  influence  of  elec¬ 
tricity.  The  zero  here  signifies  the  gas  concentration  would  not 
change  and  no  heterogeneous  chemical  reactions  would  occur. 
When  eCa  or  dca  increases,  the  effective  diffusivity  of  gas  species  Deff 
will  increase  and  the  mass  transport  process  is  improved  according 
to  Eqs.  (20)— (23),  although  the  original  mass  transfer  rate  at  the 
cathode  surface  decreases.  eca  could  affect  both  the  binary  molec¬ 
ular  diffusion  and  Knudsen  diffusion  coefficient  while  dca  could 
only  affect  the  latter  coefficient,  and  then  eca  and  dca  both 
influence  the  molecular  diffusion  coefficient  and  further  the  gas 
diffusivity  through  changing  the  gas  molar  fraction  x.  That's  why 
zone  size  of  heterogeneous  chemical  reactions  increases  with  eca 
and  dca,  and  the  variation  of  zone  size  caused  by  eca  is  bigger  than 
that  caused  by  dca,  when  eca  and  dca  are  both  doubled. 

In  addition,  the  results  in  Fig.  6(b)  and  (e)  show  that  the  charge 
transfer  flux  of^VVjon  also  drops  down  as  expected  when  away 
from  electrolyte.  The  rate  of  descent  for  charge  transfer  is  slower 
than  that  for  mass  transfer,  so  that  the  main  zone  of  electro¬ 
chemistry  is  larger  than  that  of  chemistry.  As  eca  or  dca  increases, 
the  TPB  active  area  per  unit  volume  Sj pB  decreases  and  charge 
transfer  is  affected  according  to  Eq.  (15),  which  slightly  enlarges 
the  zone  size  of  electrochemistry  like  the  small  variation  in 
Fig.  6(f).  Moreover,  the  porosity  is  related  to  the  effective  ionic 
conductivity  creff: 

tf'eff  =  (1  -  e)<r  (28) 

The  impact  of  ionic  conductivity  would  be  larger  than  that  of 
reaction  active  sites,  so  that  charge  transfer  flux  obviously  de¬ 
creases  with  £ca  and  electrochemistry  reacts  in  a  smaller  reaction 
zone. 

On  the  whole,  as  eca  or  dca  increases,  the  mass  transport  is 
obviously  improved  and  the  main  zone  of  heterogeneous  chemical 
reactions  is  expanded,  so  that  more  CO2  reacts  through  chemical 
reaction  pathway  for  CO2/H2O  co-electrolysis  and  then  the  relative 
cell  performance  is  prone  to  H2O  electrolysis. 

4.2.  Effect  of  ionic  conductivity 

The  cathode  ionic  conductivity  <7ion,ca  is  another  typical  mate¬ 
rial  parameters.  Fig.  7  compares  the  relative  polarization  curves  of 
H2O  electrolysis,  CO2  electrolysis  and  CO2/H2O  co-electrolysis  and 
Fig.  8  shows  the  heterogeneous  chemical/electrochemical  reaction 
main  zones,  with  different  multiplying  coefficients  of  ionic  con¬ 
ductivity  at  the  same  ratio  of  oxidant  and  reductant.  The  absolute 
cell  performances  of  the  three  electrolysis  modes  are  all  enhanced 


by  improving  the  ionic  conduction  or  charge  transport  as  excep¬ 
ted.  It  can  be  observed  that  the  ionic  conductivity  also  plays  an 
important  role  on  the  relative  performance  of  the  three  electrol¬ 
ysis  modes.  When  <7ion,ca  increases,  the  curve  of  co-electrolysis 
gradually  approaches  to  that  of  CO2  electrolysis,  and  more  CO2 
participates  in  electrolysis  as  shown  in  Table  4  rer/c r  at  1.4  V  in¬ 
creases  from  2.29%  to  2.60%  when  ionic  conductivity  grows  to  4 
times  larger.  This  could  be  because  that  when  the  electrochemical 
reaction  rate  increases,  the  competitive  ability  of  electrochemistry 
is  relatively  improved,  so  that  more  CO2  reacts  through  the  elec¬ 
trochemical  pathway. 

Besides,  (7ion>ca  almost  has  no  effect  on  mass  transfer  flux  at 
OCV  as  shown  in  Fig.  8(a),  so  that  the  zone  size  of  heterogeneous 
chemical  reaction  doesn't  change  with  <7ion,ca  in  Fig.  8(c).  The 
results  also  show  that  <7ion,ca  significantly  increases  the  charge 
transfer  flux  and  expands  the  electrochemical  reaction  main 
zone.  As  (7ion,ca  increases  to  4  times,  the  size  of  electrochemical 
main  zone  is  enlarged  from  438  pm  to  477  pm.  Therefore,  the 
results  above  indicate  that  when  the  zone  of  electrochemistry  is 
enlarged,  more  CO2  for  CO2/H2O  co-electrolysis  could  be  elec¬ 
trolyzed,  resulting  in  a  polarization  curve  approaching  to  the 
curve  of  CO2  electrolysis. 

43.  Effect  of  gas  composition 

If  only  the  ratios  of  oxidant  and  reductant  are  the  same,  the 
polarization  curves  for  the  three  electrolysis  modes  could  be 
comparable.  At  a  fixed  ratio  of  oxidant  and  reductant,  increasing 
H2O/CO2  makes  the  performance  close  to  H2O  electrolysis  as  ex¬ 
pected.  Thus  the  H2O/CO2  is  kept  as  1  and  only  the  different  oxi¬ 
dant:  reductants  are  discussed  in  details  in  this  section.  The  relative 
polarization  curves  and  heterogeneous  chemical/electrochemical 
main  zones  are  displayed  in  Figs.  9  and  10. 

The  results  in  Fig.  9  and  Table  4  show  that  when  oxidant: 
reductant  =  4:1,  the  polarization  curves  of  CO2/H2O  co-electrolysis 
and  H2O  electrolysis  are  very  close,  and  only  0.35%  CO2  could  be 
electrolyzed  at  1.4  V.  As  oxidant:  reductant  increases,  the  relative 
polarization  curve  of  C02/H20  co-electrolysis  becomes  prone  to 
that  of  CO2  electrolysis.  Meanwhile,  the  ratio  of  CO2  participating  in 
electrolysis  and  heterogeneous  reactions  rer/c r  at  1.4  V  significantly 
grows  from  0.35%  to  12.03%.  The  tendency  of  relative  polarization 
curve  is  in  accordance  with  the  variation  of  rer/c r.  This  variation  of 
rer/cr  is  much  bigger  than  that  caused  by  material  characteristics  of 
cathode  in  the  range  of  this  simulation. 

Because  the  diffusion  coefficient  of  H2  is  higher  than  the  other 
3  gas  species,  bigger  H2  molar  fraction  in  gas  composition  should 


Fig.  9.  Polarization  curves  comparison  of  H20  electrolysis,  C02  electrolysis  and  C02/H20  co-electrolysis  with  different  ratios  of  oxidant  and  reductant. 
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contribute  to  higher  effective  diffusivity  and  better  mass  transfer 
without  regard  to  reactions.  However,  from  Fig.  10(c),  the  zone 
size  of  heterogeneous  chemical  reactions  doesn't  behavior  a 
linear  variation.  The  results  in  fable  5  show  that  the  original 
effective  diffusivities  of  the  4  species  all  increase  with  the  H2 
molar  fraction  at  the  cathode  surface  x  =  0  pm  at  OCV,  as  the 
speculation  above.  However,  at  the  interface  of  electrolyte  and 
cathode  x  =  695  pm,  the  diffusivities  of  H2  decrease  first  and  then 
increase  while  that  of  CO2  increase  first  and  then  decrease.  This 
phenomenon  is  mainly  attributed  to  the  heterogeneous  chemical 
reactions,  which  highly  control  the  variation  of  gas  concentration 
Vc  and  then  further  influence  Deff.  This  suggests  that  to  analyze 
the  mass  transfer,  the  diffusivities  Deff  should  be  taken  into 
consideration,  as  well  as  the  variation  of  gas  concentration  Vc 
caused  by  reactions.  According  to  the  mass  transfer  flux  in 
Fig.  10(a),  the  heterogeneous  chemical  reactions  make  D^2VcCq2 
at  oxidant: reductant  =  4:1  approaches  to  zero  with  the  fastest 
rate  and  that  at  oxidant:reductant  =  2:3  varies  to  zero  the  most 
slowly.  These  results  mean  that  the  mass  transfer  is  improved 
when  oxidant:  reductant  varies  from  4:1  and  2:3,  and  then  is 
reduced  when  oxidant:  reductant  varies  from  2:3  to  1:4,  with 
both  diffusivities  and  reactions  taken  into  account.  Therefore,  the 
zone  size  of  heterogeneous  chemical  reactions  is  enlarged  first 
and  then  is  reduced  as  shown  in  Fig.  10(c). 

On  the  other  hand,  although  the  absolute  electrochemical  re¬ 
action  rate  and  the  original  charge  transfer  flux  at  x  =  695  pm  are 
reduced  with  oxidant:  reductant,  the  charge  transfer  fluxes 
of^Wjon  of  the  5  cases  show  a  totally  opposite  behavior  when  x  is 
close  to  0  as  shown  in  Fig.  10(b):  of0fnVVion  gradually  increases  with 
oxidant:  reductant.  Thus,  the  zone  size  of  electrochemical  reactions 
significantly  increases  from  319  pm  to  549  pm.  Overall,  the  gas 
composition  both  significantly  affects  the  heterogeneous  chemical/ 
electrochemical  reaction  main  zones.  The  influence  on  electro¬ 
chemistry  is  bigger  than  that  on  heterogeneous  chemistry,  so  that 
the  overall  relative  cell  performance  is  prone  to  C02  electrolysis 
when  oxidant:  reductant  increases. 


0  100  200  300  400  500  600  695 

Distance  from  outside  surface  of  cathode  (pm) 

Fig.  10.  Effects  of  the  ratio  of  oxidant  and  reductant  on  heterogeneous  chemical / 
electrochemical  reaction  main  zones  in  cathode  for  C02/H20  co-electrolysis. 


4.4.  Effect  of  temperature 

As  another  important  operating  parameter,  the  effect  of  tem¬ 
perature  on  the  relative  polarization  curves  of  H2O  electrolysis,  CO2 
electrolysis  and  C02/H20  co-electrolysis  and  on  the  heterogeneous 
chemical/electrochemical  reaction  main  zones  are  displayed  in 
Figs.  11  and  12.  Temperature  could  significantly  increase  the  ab¬ 
solute  cell  performance  of  the  three  electrolysis  modes,  and  also 
could  make  the  polarization  curve  of  CO2/H2O  co-electrolysis 
approach  to  that  of  CO2  electrolysis.  It  should  be  noted  that  the 
polarization  curve  of  CO2/H2O  co-electrolysis  in  Fig.  11(c)  is  worse 
than  that  of  CO2  electrolysis  at  lower  voltage,  just  because  the  OCV 
of  co-electrolysis  at  900  °C  is  higher  than  that  of  C02  electrolysis.  At 
higher  voltage,  the  curve  of  co-electrolysis  lies  between  the  other 
two  electrolysis  modes.  The  results  in  Table  4  also  represent  that 


Table  5 

Effective  diffusivity  of  gas  species  with  different  ratios  of  oxidant  and  reductant  for  C02/H20  co-electrolysis  at  OCV. 


Oxidant:  reductant 

<o(10- 

4  m2  s"1) 

(10-4  m2 

s'1) 

D^2  (10-4  m2  s-1) 

Dg  (10-4  m2  s"1) 

x  =  0  pm 

x  =  695  pm 

x  =  0  pm 

x  =  695  pm 

x  =  0  pm 

x  =  695  pm 

x  =  0  pm 

x  =  695  pm 

4:1 

1.03 

1.03 

2.13 

2.10 

0.77 

0.83 

0.863 

0.85 

3:2 

1.10 

1.10 

2.14 

2.04 

0.79 

0.94 

0.913 

0.89 

1:1 

1.14 

1.13 

2.14 

2.01 

0.80 

0.99 

0.946 

0.91 

2:3 

1.17 

1.16 

2.15 

2.01 

0.81 

0.98 

0.987 

0.94 

1:4 

1.23 

1.21 

2.15 

2.08 

1.23 

0.87 

1.107 

1.04 
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Fig.  11.  Polarization  curves  comparison  of  H20  electrolysis,  C02  electrolysis  and  C02/H20  co-electrolysis  with  different  temperatures. 


more  CO2  could  be  electrolyzed  when  increasing  temperature.  rer/c r 
at  1.4  V  significantly  grows  from  2.29%  to  8.23%. 

As  temperature  rises  up,  diffusion,  ionic  conduction,  rates  of 
chemical  and  electrochemical  reactions  could  be  all  enhanced.  So 
according  to  Fig.  12(a)  and  (b),  although  the  original  Vcco2 
x  =  0  pm  decreases  with  temperature,  both  mass  and  charge 
transfer  are  improved  by  temperature,  especially  for  charge 
transfer.  The  sizes  of  heterogeneous  chemical/electrochemical 
reaction  main  zones  are  both  enlarged  with  temperature,  as 
shown  in  Fig.  12(c).  Besides,  the  electrochemical  reactions  are 
more  sensitive  to  temperature  than  chemical  reactions,  so  that  as 
temperature  rises  up,  relatively  more  CO2  reacts  through  elec¬ 
trolysis  and  then  the  polarization  curve  of  co-electrolysis  is  prone 
to  CO2  electrolysis. 

5.  Conclusions 

In  this  paper,  a  one-dimensional  elementary  reaction  model  of 
solid  oxide  electrolysis  cells  coupling  with  the  heterogeneous 
elementary  reactions,  electrochemical  reactions  and  the  mass  and 
charge  transport  is  applied  to  study  the  relative  performance  of  H20 
electrolysis,  CO2  electrolysis  and  CO2/H2O  co-electrolysis  and  the 
competitive  behavior  of  heterogeneous  chemical  and  electro¬ 
chemical  reactions.  The  two  kinds  of  reactions  exist  in  different  zones 
in  the  cathode  of  SOECs:  the  heterogeneous  chemical  reactions 
mainly  react  near  the  outside  surface  of  cathode,  while  the  electro¬ 
chemical  reactions  mainly  occur  near  the  electrolyte.  The  sizes  of  the 
heterogeneous  chemical/electrochemical  reaction  main  zones  are 
quantified  according  to  the  mathematical  analysis.  The  effects  of 
material  characteristics  of  cathode  (porosity,  particle  diameter  and 
ionic  conductivity)  and  operating  conditions  of  SOECs  (gas  compo¬ 
sition  and  temperature)  are  studied  in  details.  All  the  comparisons  of 
the  relative  cell  performance  require  that  the  ratios  of  oxidant  and 
reductant  of  the  three  electrolysis  modes  should  be  the  same. 

On  the  whole,  the  modeling  results  indicate  that: 

(1)  The  rate  of  heterogeneous  chemical  reactions  is  about 
20-100  times  faster  than  that  of  electrochemical  reactions. 
At  the  same  local  reaction  site,  if  the  reactant  gas  doesn't 
reach  the  equilibrium,  the  gas  prefers  to  heterogeneous 
chemical  reactions  rather  than  electrochemical  reactions. 

(2)  The  mass  transfer  flux  DVc  determines  the  main  zone  size  of 
heterogeneous  chemical  reactions,  while  the  charge  transfer 
flux  (jW  determines  the  main  zone  size  of  electrochemical 
reactions. 


(3)  When  the  mass  transfer  is  improved,  the  zone  size  of  het¬ 
erogeneous  chemistry  could  be  enlarged,  more  CO2  could 
react  through  heterogeneous  chemical  pathway,  and  polari¬ 
zation  curves  of  CO2/H2O  co-electrolysis  could  be  prone  to 
that  of  H2O  electrolysis. 

(4)  When  the  charge  transfer  is  improved,  the  zone  size  of 
electrochemistry  could  be  enlarged,  more  CO2  could  react 
through  electrochemical  pathway,  and  polarization  curves  of 
CO2/H2O  co-electrolysis  could  be  prone  to  that  of  CO2 
electrolysis. 

Porosity,  particle  diameter,  ionic  conductivity,  gas  composition 
and  temperature  have  different  influence  on  mass  and  charge 
transfer  in  cathode,  and  then  further  differently  affect  the  hetero¬ 
geneous  chemical/electrochemical  reaction  main  zones  and  the 
relative  cell  performance  of  the  three  electrolysis  modes. 

(1)  Porosity  and  particle  diameter  strongly  affect  the  mass 
transfer  flux  and  the  zone  size  of  heterogeneous  chemical 
reaction,  but  slightly  affects  the  charge  transfer  flux  and  the 
zone  size  of  electrochemical  reaction.  As  porosity  or  particle 
diameter  increases,  the  mass  transport  is  obviously 
improved  and  the  main  zone  of  heterogeneous  chemical 
reactions  is  expanded,  so  that  more  CO2  reacts  through 
chemical  reaction  pathway  for  CO2/H2O  co-electrolysis  and 
then  the  relative  cell  performance  of  co-electrolysis  is  prone 
to  that  of  H2O  electrolysis. 

(2)  Increasing  ionic  conductivity  has  no  effect  on  mass  transfer 
flux,  but  significantly  improves  charge  transfer  flux  and  en¬ 
larges  the  zone  size  of  electrochemical  reaction.  Thus,  the 
competitive  ability  of  electrochemistry  is  relatively 
improved,  more  CO2  could  be  electrolyzed,  and  polarization 
curve  of  CO2/H2O  co-electrolysis  approaches  to  the  curve  of 
C02  electrolysis. 

(3)  Gas  composition  shows  the  biggest  influence  on  the  two 
zones  with  the  parameters  in  the  range  of  this  simulation. 
The  zone  size  of  heterogeneous  chemical  reactions  is 
enlarged  first  when  oxidant:  reductant  increases  from  4:1  to 
2:3,  and  then  is  reduced  when  oxidant:  reductant  varies  from 
2:3  to  1:4.  The  zone  size  of  electrochemical  reactions 
significantly  increases  with  oxidant:  reductant.  The  influence 
on  electrochemistry  is  bigger  than  that  on  heterogeneous 
chemistry,  so  that  when  oxidant:  reductant  increases,  quite 
more  CO2  could  react  through  electrochemistry  and  the 
overall  relative  cell  performance  of  co-electrolysis  is  prone  to 
CO2  electrolysis. 
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Fig.  12.  Effect  of  temperature  on  heterogeneous  chemical/electrochemical  reaction 
main  zones  in  cathode  for  C02/H20  co-electrolysis. 

(4)  As  temperature  rises  up,  both  mass  and  charge  transfer  are 
improved  and  the  sizes  of  two  zones  expand.  The  electro¬ 
chemical  reactions  are  more  sensitive  to  temperature  than 
chemical  reactions,  thus  when  temperature  increases,  rela¬ 
tively  more  CO2  reacts  through  electrolysis  and  the  polariza¬ 
tion  curve  of  co-electrolysis  is  prone  to  that  of  CO2  electrolysis. 
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Nomenclature 


Abbreviation 


LSM:  lanthanum  strontium  manganite 
RWGS:  reverse  water  gas  shift 
ScSZ:  scandium  stabilized  zirconium 
SOEC:  solid  oxide  electrolysis  cell 
TPB:  triple  phase  boundary 
WGS:  water  gas  shift 
YSZ:  yttrium  stabilized  zirconium 

English  letter 


a:  pre-exponential  factor  in  sticking  coefficient  expression 
A:  pre-exponential  factor  of  the  Arrhenius  form  (in  terms  of  cm,  mol  and  s) 
b:  temperature  exponent  in  sticking  coefficient  expression 
c;  concentration  of  species  (mol  m-3) 

C(Niy.  surface  concentrations  of  the  free  surface  active  sites  on  the  Ni  surface 
(mol  m-2) 

co(Nif  surface  concentrations  of  O  element  on  the  Ni  surface  (mol  m-2) 
co2  (ysz)  :  the  volumetric  concentrations  of  interstitial  oxygen  in  the  YSZ  ionic 
conductor  (mol  m-3) 

C(ysz)'  the  volumetric  concentrations  of  interstitial  oxygen  in  the  YSZ  ionic  conductor 
(mol  m-3) 

CquII<:  the  oxygen  concentrations  in  the  bulk  (mol  m-3) 
c(f2B:  the  oxygen  concentrations  at  the  TPB  (mol  m-2) 
d;  activation  energy  in  sticking  coefficient  expression  (J  mol-1) 

Deff;  effective  diffusion  coefficient  (m2  s-1) 

D^:  effective  Knudsen  diffusion  coefficient  (m2  s-1) 

Dynoie'  effective  molecular  diffusion  coefficient  (m2  s-1) 

E;  activation  energy  (kj  mol-1) 

F:  Faraday  constant  (96,384C  mol-1) 
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io:  exchange  current  density  (A  m-2) 
k:  reaction  rate  constant  (in  terms  of  m,  mol  and  s) 
kec:  forward  electrochemical  reaction  rate  (mol  m-2  s-1) 
k-ec:  reverse  electrochemical  reaction  rate  (mol  m-2  s-1) 

Kg:  number  of  gas  species 

I<s:  number  of  surface  species 

M:  molecular  weight  (kg  mol-1) 

n;  reaction  order  of  the  Arrhenius  form 

nt:  total  number  of  electronic  and  ionic  conductor  particles 

nep:  fraction  number  of  electronic  conductor  particles 

nip:  fraction  number  of  ionic  conductor  particles 

N:  number  of  reactions 

p;  pressure  (Pa) 

Pep:  whole  range  connection  probabilities  of  electronic  conductor  particles 
Pip-  whole  range  connection  probabilities  of  ionic  conductor  particles 
r:  average  pore  radius  (m) 

rep:  mean  radius  of  the  electronic  conductor  particle  (m) 

R:  gas  constant  (8.314  mol-1  K-1)  or  source  term  of  mass  balance  equation 
(kg  m-3  s-1) 

Se effective  reaction  area  per  unit  volume  (m2  m-3) 

Sm-  Ni  active  surface  area  per  unit  volume  (m2  m-3) 

Stpb :  TPB  active  area  per  unit  volume  (m2  m-3) 

T:  temperature  (K) 

Q;  current  source  (A  m-3) 

V:  diffusion  volume  or  electronic  potential  (V) 

W:  molecular  weight  of  gas  species  (kg  mol-1) 
x:  molar  fraction 

Z;  mean  coordination  number  of  electron  and  ionic  conductor  particles 
Zep:  coordination  number  of  electron  conductor  particles 
ZiP:  coordination  number  of  ionic  conductor  particles 

Greek  letters 

a:  charge  transfer  coefficient 

(3:  cathode  electrochemical  kinetics  parameter 

y:  the  sum  of  all  of  the  surface  reactants  stoichiometric  coefficients 

T:  surface  sites  density  (Q-1  m-2) 

e :  porosity  or  parameter  modeling  the  species  coverage 

7):  overpotential  (V) 


0:  surface  coverage  or  contact  angle  between  the  electronic  and  ionic  conductors 
(rad) 

p:  parameter  modeling  the  species  coverage 
v':  stoichiometric  coefficient  of  the  reactants 
v  :  stoichiometric  coefficient  of  the  products 
aeff;  corresponding  conductor  phase  effective  conductivity  (S  m-1) 
r:  tortuosity 
species  symbol 
dQ:  computational  domain 

Subscripts 

ac:  anode  chamber 
act:  active  layer 
an:  anode 
ca:  cathode 

cccr:  cathode  chamber  chemical  reactions 
ec:  electrochemical  reactions 
el:  electrolyte 
elec:  electronic 

ep:  electronic  conductor  particle 
er:  electrochemical  reactions 
i:  reactions  index 
g:  gas  phase  species 
ion:  ionic 

ip:  ionic  conductor  particle 
Kn:  Knudsen 
mole:  molecular 
ref:  reference 
sp:  support  layer 
t:  total 

Superscripts 

0:  parameter  at  equilibrium  conditions 
bulk:  bulk  phase 
eff:  effective 

TPB:  triple  phase  boundary 


